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We discuss of a ring-shaped soft quantum wire modeled by 5 inter- 
action supported by the ring of a generally nonconstant coupling 
strength. We derive condition which determines the discrete spec- 
trum of such systems, and analyze the dependence of eigenvalues and 
eigenfunctions on the coupling and ring geometry. In particular, we 
illustrate that a random component in the coupling leads to a localiza- 
tion. The discrete spectrum is investigated also in the situation when 
the ring is placed into a homogeneous magnetic field or threaded by 
an Aharonov-Bohm flux and the system exhibits persistent currents. 



1 Introduction 



The aim of the present paper is to investigate spectral properties of a two- 
dimensional quantum particle subject to a 5 interaction supported by a cir- 
cular ring, plus possibly a magnetic field. Similar systems were studied for 
a quite a long t ime. On e can mentioned an early paper by Kurilev |Kuj and 
later work jRTl iREKSj devoted to measure-type perturbations of the Lapla- 
cian where some general spectral results were derived 1 . Later this line of 

1 Recall also a related model of a photonic crystal studied in [FKUKuK] . It differs from 
the present by the switched roles of the coupling and spectral parameters. 
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study was extended to the situation when a singular zero-measure set itself 
supports a nontrivial dynamics \K&\ IKK| IKoj . 

In the most part of the mentioned work the emphasis was laid on methods 
to treat such singular interactions. Recently the problem was put into a 
different perspective related to a specific physical model describing idealized 
"quantum wires" . These tube-like or graph-like semiconductor structures are 
usually modeled as strips or tubes with hard walls, or in a more simplified 
version by graphs supporting "one- dimensional" electrons. A weak point 
of such models is the assumption about the nature of the confinement. In 
actual quantum wires, the electrons are trapped due to interfaces between two 
different semiconductor materials which represents a finite potential jump. 
Thus electrons can be found also outside such a wire, although not "too far" 
since the exterior represents here a classically forbidden region. 

Recent investigations revealed interesting relations between the spec- 
trum and geometry of such a "soft" quantum wire, namely the existence of 
curvature- induced bound states |Elj . Moreover, it was demonstrated that in 
the strong-coupling limit the negative part of the spectrum of such Hamilto- 
nians approaches that of an "ideal" curve-like wire with an effective potential 
determined by the curvature |EYH lEY2j . On the other hand, apart of a few 
very simple examples we are lacking solvable models of such systems. 

In the present paper we analyze a simple model with a S interaction 
supported by a circle of radius R and a generally non-constant coupling 
strength. In particular, we are interested in the discrete spectrum of such ring 
systems. First we illustrate that it is rather the geometry of the interaction 
curve than its topology which determines the spectral properties. To this aim 
we discuss in Section 2 the case of the "full" and "broken" ring, the latter 
being obtained by putting the coupling constant zero at its segment; we will 
show that the spectrum depends substantially on the ability of the particle 
to tunnel though the "gap". Another natural question concerns spectral 
properties in the situation when the coupling strength is randomly varying. 
It is customary to investigate it in an infinite-length setting as a passage from 
absolutely continuous to dense pure point spectrum [St]. At the same time 
the localization effect can be seen in systems with perturbations of a finite 
length too if we observe the characteristic size of the wavefunction; we will 
illustrate this fact in the framework of the present model. 

Another situation we are going to discuss concerns the case when the 
particle is exposed to a magnetic field perpendicular to the ring plane. In 
Section 3 we discuss two models, one with a homogeneous background field 
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and the other with a pointlike Aharonov-Bohm flux piercing the centre of the 
ring. For an ideal wire these situations are, of course, equivalent because the 
only quantity which counts is the flux through the loop. This is no longer 
true in the soft case when the particle wave packet extends outside the ring. 

An important question concerns the existence of current-carrying states. 
Recall that persistent currents in rings threaded by a magnetic flux are one 
of the characteristic features of mesoscopic systems - see, e.g., fCGKl fCWB 
and numerous other theoretical and experimental papers where they were 
discussed. If an electron is strictly confined to a loop T the effect is con- 
tained in the dependence of the corresponding eigenvalues Ej on the flux 
threading the loop (measured in the units of flux quanta, 2nhc\e\~ 1 ); the 
persistent current Ij in the j-th state is defined as the multiple —cdEj/dcf) of 
the eigenvalue derivative w.r.t. the flux. In particular, if the particle motion 
on such a loop is free, we have 

where L is the loop circumference, so the currents depend linearly on the 
applied field. For a soft quantum wire in the form of a closed curve with an 
attractive coupling which is constant around the loop and strong enough the 
existence of persistent currents was demonstrated in |EY3j . An illustration 
of this effect can be seen in a related model |CE] in which a curve supported 
(5-coupling is replaced by an array of two-dimensional point interaction. In 
analogy with this paper we also expect that a randomization of the coupling 
will destroy a coherent transport around the ring, however, since this problem 
is numerically very demanding we postpone it to a subsequent paper. 

2 Formulation of the problem 

Let r := {x : \x\ = R} with an R > be a circle in M. 2 and a : r — > R a 
piecewise continuous function. We are going to discuss the Hamiltonian in 
L 2 (R 2 ) which can be formally written as 

tfform = HV - A(r)) 2 - a(<p)5(r - R) , (2.1) 

where r, ip are the polar coordinates of a point x and A is a vector potential 
determining the magnetic field. We suppose that the latter is rotationally 



3 



symmetric. For convenience we associate positive values of the coupling 
strength a with the attractive interaction. As usual we also simplify the 
treatment by using rationalized units, h = 2m* = e = c = 1. 

The operator (|2.1j) can be defined properly, and quite generally, by means 
of the corresponding quadratic form. As long as the function a is suffi- 
ciently regular, however, one use alternatively boundary conditions BEKSj 
Rem. 4.1], i.e. to consider the operator H a ^ acting as 



(H a , A ip) (x) = ((-iv - A) 2 i[>) (x) , x e m 2 \ r , 

for any ip of the domain consisting of functions which belong to 1^2,2(1 
are continuous at V with the normal derivatives having a jump there. 



9fp 
dr 



r=R+ 



dip _ 
dr 



-a((p)ip(x) , x = (r, Lp) E T . 



(2.2) 

e\r), 

(2.3) 



It is straightforward to check that H 0) a is e.s.a., so the Hamiltonian of our 
system corresponding to the formal expression (|2.1jl can be identified with 
its closure. 



2.1 A warm-up: constant 5 interaction on a ring 

Let us consider first the simplest case when the magnetic field is absent, the 
coupling strength is constant along the ring and the interaction is attractive, 



a(cp) = a > 0. We look for the discrete spectrum which is in view of [BEKSl 
Thm. 4.2] located at the negative halfline. Since the system is rotationally 
symmetric we can perform the partial-wave decomposition and adopt the 
following Ansatz for the eigenfunctions 

ip m (r^)=Pm(r)e vmip } meZ. 

Outside the ring the solution coincides with that of the free Schrodinger 
equation, including the behaviour at the origin and at infinity, i.e. 



Pm{r) 



Ci/ TO («r) ... r < R 
c 2 K m (Kr) ... r > R 



corresponding to the energy — k 2 . Using the continuity of the function p m at 
r = R together with the matching condition (|2.3|) and the Wronskian relation 
W(K m (z), I m ( z )) = w e get the spectral condition 

{I m K m ){KR) = — . (2.4) 
aR 
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To find its solutions, recall some properties of the modified Bessel functions 
|ASt Chap. 9]; since I_ m (z) = Im(z) an d K-m{z) — K m (z), we may consider 
m > only. The l.h.s. of (|2.4j) is decreasing in (0, oo) and has the following 
asymptotics, 



-In | -7 + 0(1) ... m = 
^ + 0(z) ... m^O 

as z — > 0, where 7 = 0.577 ... is the Euler number, and 



2 K-i) lC v^ 



as z — > 00. Consequently, there is at most one solution at each partial wave 
which we have anticipated labeling the solution by the orbital index only. 
The ground state corresponding to m = exists for any a > 0, while other 
partial waves exhibit a bound state provided 

aR > 2\m\ . 

In the weak-coupling case, aR -C 1, the binding energy is exponentially 
small, 

as expected for a two-dimensional Schrodinger operator 2 , while for aR 3> 1 
the above asymptotics gives 

n 2 ™2 — 1 

E m = ~\ + ^ + 0(a-*R-*) . 

This can be regarded as a particular case of the theorem proven in |EYlj . just 
the error estimate is better than in the general case. The discrete spectrum 
for a fixed coupling strength as a function of R is plotted in Fig. 1. 

One can continue the analysis of the rotationally symmetric (^-interaction 
on the ring by discussing the scattering, approximations by means of scaled 
potentials, combinations with other interactions, etc., in analogy in the three- 
dimensional case considered in |AGS| . We will not do that, however, because 
we are interested mainly in the situation without the rotational symmetry. 

2 Taking the limit R — > with a properly scaled a{R) we can get a two-dimensional 
point inter action. If a depends on the angle, such a limit can be done using the resolvent 



formula of BEKSJ in analogy with the tree-dimensional case worked out recently in |Sh| 
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2.2 General formalism 



Next we want to look what happens if the rotational symmetry is broken 
due to non-constancy of the function cn(-). In this case the point interaction 
couples different partial waves and we have to look for the solution in the 
form of a series. Since the scheme is similar for the magnetic case, or even 
for operators 

H aA + V(r) 

with a rotationally invariant potential, we shall formulate it generally adopt- 
ing the following Ansatz 

where {c m }, {d m } are coefficient sequences and f m , g m are solutions of the 
"free" Schrodinger equation with the energy — k? inside and outside the ring, 
respectively, in the m-th partial wave; recall that by assumption the system 
without the point interaction is rotationally symmetric. We impose the stan- 
dard requirement that the solution is regular at the origin and L 2 at infinity, 
then the functions f m , g m are unique up to multiplicative constants which 
can be absorbed into the coefficients. 

The above function has naturally to belong to the domain of the Hamil- 
tonian. Its continuity at r = R together with the orthonormality of the 
trigonometric basis in L 2 ((0,27r)) implies 

Cm = d m ——— . (2.6) 

Jm{rC) 

At the same time, the matching condition (|2.Hj) yields 



imip 



multiplying this relation by e~ mif and integrating w.r.t. ip we get an ex- 
pression of c n by means of the coefficients d m . Combing it with ()2.6|) and 
introducing 



a 



mn •" 



^ a(<p) e^ m ~ n ^ dp , (a) := a mm = [** a(<p) dp 
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we get after a simple manipulation 

5^ {Smn \(o)fn9n + 2irW (f n , g n )\ (R) + (l - 8 mn ) a mn (g m f n ) (R) \d m = 

(2.7) 

as a set of equations for the coefficients, or in other words, an operator 
equation in the £ 2 space of the coefficients. Solving it numerically, one takes 
a family of truncated systems; the con vergence of such an approximation is 
checked in the same way as in [E STVj . 



2.3 A broken ring 

Consider now the case when the coupling strength is a step function, 

®(<P) = <XX16/2,2k-$/2]('P) 

for some a > and 9 G (0, 2n). Substituting W(I m (KR), K m {nR)) = -R~ l 
(recall that the derivative is taken w.r.t. r) and 



a, 



a(2n — 9) ... m = n 

2a • 0(m-n) / 

sin -H; — - ... m t 1 n 

m-n 2 ' 



into (|2.7j) we get equations for the coefficients d m in this case. 

In an ideal quantum wire described by an appropriate Schrodinger op- 
erator on T it is the topology rather than geometry which determines the 
character of the discrete spectrum. For the full ring, 9 = 0, the spectrum is 
twice degenerate with the exception of the ground state, while a broken ring 
spectrum looks like that of an infinite deep rectangular well. The spectrum 
of a soft wire will be similar if the distance between the loose ends is large 
so that the tunneling between them is negligible. To get an idea at what 
values of the parameters the regime changes, recall that the eigenfunction 
of the one-center point interaction of the coupling strength —a on line is 
e -«M/ 2 [AGHHj . and thus the characteristic "size" of the bound-state wave 
function measured from the ring is 2ct -1 . Comparing this to the "gap size" 
we see that the switch from the full-ring-type spectrum to the infinite-well- 
type takes place at a ~ 29R. This is illustrated in Fig. 2 where we plot the 
spectrum as a function of the gap angle; it is seen that the transition takes 
place around the value where we expect it. 
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3 Magnetic rings 

3.1 The homogeneous background field 

Consider the magnetic field of intensity B > perpendicular to the ring 
plane. In the circular gauge the corresponding vector potential is A — | Br e v 
and the Hamiltonian without the singular interaction at the ring acts as 

d 2 Id 1 / . d 1 „ 2 N 2 



After a standard manipulation we find that the solutions at energy z for an 
interior and exterior of the ring, to be inserted into (|2.5|) . are 

f m (r) = rH e -^ 2 /4 M (a m (z),\m\ +l;^Br 2 ^J , 



9m(r) 

respectively, where 



I \ 1 ( 2 



and M, U are the regular and singular confluent hypergeometric function, 
respectively. In particular, the spectrum of the free (Landau) Hamiltonian 
is given by the condition a m (z) = —n, n = 0, 1,2, . . ., i.e. consists of the 
Landau levels 

z = B(m + \m\ + 2n + 1) , n = 0,1,2,.... 

To employ f)2.7j) we need also the Wronskian of the above solutions which is 
by HSl 13.1.22] equal to 



r T(a m (z)) \2 

In particular, if the ring is rotationally symmetric, oe(ip) = a, the condition 
(|2.7|) is reduced to 

(1 p>f>2\\ m \ / i \ 1 

9l II na m {z))e~ BR2 /\MU) la m (z), H + 1;~BR 2 ) = — - . 
2\m\\ \ 2 J an 
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3.2 The AB flux case 

Suppose now that instead of the homogeneous field considered above the ring 
is threaded at its centre by a magnetic flux line. The vector potential can be 
now chosen as A = e„, where $ is the value of the flux. The normalized 
flux appearing in (Jl.l)) equals in the present units <fi = $/2-7r. 

We a ssume that the flux line is not combined with a point interaction 
|AT| IDS] ; treating then the Schrodinger equation in a standard way [Ruj we 
get the interior and exterior solutions at the energy — k, 2 in the form 

f m (r) = I\ m ^(Kr) , g m (r) = K\ m ^\(Kr) , 

respectively; their Wronskian is 

1 



W(f m (r),g m (r)) 



r 



In particular, in case of a rotationally symmetric ring we get the spectral 
condition which is similar to (|2.4)1 . namely 

(J| m _0|iC| m _^l)(«i2) = ^= . (3.1) 
Taking into account the relation 

1 

\im(I\ u \K\ u \)(z 



we can find the condition 

2\m-(j)\ = aR (3.2) 

that determines critical values at which the eigenvalues are absorbed in the 
essential spectrum. 



4 Results and discussion 

Let us now turn to a more detailed discussion of the results. 
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4.1 Non- magnetic case, constant a 



We have presented already in Figs. 1,2 the spectrum for the full and broken 
ring, respectively. If the radius is large enough (in terms of the quantity aR) 
the spectrum has a distinctive one-dimensional character with the largest 
density at its bottom. At the same time, a gap which is large enough to 
prevent tunneling between the loose ends gives the spectrum the hard-wall 
nature. This leads us to the conjecture that the asymptotic formula derived in 
|EYlj will hold for any finite, smooth, and non-closed curve as well, with the 
comparison operator S = —-£2 — \k(s) 2 , where k(s) is the curvature, being 
now specified by Dirichlet boundary conditions at the curve endpoints. The 
one-dimensional character can be seen also from the eigenfunctions: in Fig. 3 
we plot the ground state and the eighth excited state. 

4.2 Non- magnetic case, localization 

As we have said in the introduction one expects that irregular variations of 
the coupling constant may destroy the coherence of the wavefunction along 
the curve supporting the interaction. Since the negative spectrum of the 
present model is discrete in any case, one can try to see the effect through the 
characteristic size of the eigenfunctions. If the coupling strength a is constant 
over the whole ring or a substantial part of it, the latter are distributed 
roughly uniformly along T; the uniformity is ideal for the full ring and it has 
sinusoidal variations due to the natural quantization for excited states. 

When we want to find a quantity to characterize the "size" of the eigen- 
function ip^ corresponding to the energy we notice first that it can be 
expressed as jPoj 



possibly up to a normalization factor, where y(f) is a point of the ring and 
Go is the Green function of the two-dimensional Laplacian. Consequently, 
it is sufficient to characterize ipk by its restriction to ipk(<p) = ^(^(v 3 )) to 
the ring. Then we choose as the quantity of interest the second moment of 
the corresponding probability distribution, minimized over the choice of the 
reference point, 



ipk{x) = / ij(y(ip))G (x,y((p);E k )dp , 



(A-0) := min 
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Now we take a family of random step functions ot(-) and plot the above 
quantity for the ground-state eigenfunction in dependence on the dispersion 
Aq; the result is plotted in Fig. 4. It is obvious that in the average a stronger 
coupling randomness means a "more localized" eigenfunction. 

4.3 Magnetic rings 

Let us pass now to the case with a non-zero magnetic field. If the latter 
in homogeneous, the spectrum is pure point and it accumulates at the Lan- 
dau levels. In Fig. 5 we show the eigenvalues belonging to the three lowest 
"bands" of the full ring as functions of the magnetic field intensity. Fig. 6 
shows the situation in a cut ring with B fixed and the gap angle 9 changing. 
At the limit 9 — ► 2n the eigenvalues are, of course, absorbed in the Lan- 
dau levels. What is more interesting are the avoided crossings in the higher 
"bands" which witness about the non-trivial character of the cyclotronic mo- 
tion around such an "obstacle". Another illustration of the fact is given by 
the eigenfunction contour plots shown in Fig. 7. 

On the other hand, in the Aharonov-Bohm case the essential spectrum 
covers the positive halfline, and the dependence of the negative eigenvalues on 
the magnetic flux shown in Fig. 8 resembles a picture corresponding to (jl.lj) 
for the ideal ring graph, that is a family of shifted parabolas. The same can 
be said, of course, about the lower part of the homogeneous-field spectrum 
in Fig. 5, but in distinction to that case the AB spectrum is periodic modulo 
an integer number of flux quanta. A common feature is a strong dependence 
of the eigenvalues on the flux, which by the formula given in the introduction 
means that such system exhibit non-negligible persistent currents. 

Recall in this connection a close relation between such current-carrying 
states and the edge currents |Hal IMS] . Furthermore, to create a magnetic 
transport one does not need a hard wall; much more gentle perturbations 
like a potential "ditch" , translationally invariant field modification |Iwt MP , 
lEKj or even an array of point obstacles K.IK 11 are sufficient. If a "linear" 
perturbation is replaced by a circular one, one can expect occurrence of 
localized states carrying current along the circle. Naturally, the magnetic 
transport may be destroyed by a disorder. The problem of stability of (a 
part of) absolutely continuous spectrum with respect to perturbations has 
been studied recently for edge currents in halfplanes and similar domains 
[BT1 IFGWI IMMPj . strips |K,)K2| and compact domains jEMl- In a finite- 
length setting the localization due to disorder was shown in |(JKj ; we intend to 
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do the same in the context of the present model in a subsequent publication. 

While two types of the magnetic field yield similar results, there are im- 
portant differences. One of them concerns the character of the perturbation. 
The Aharonov-Bohm Hamiltonian is not an analytic perturbation of the op- 
erator with = 0. This is obvious from the fact that the eigenvalue curves 
are in this case non-smooth at the integer values of the flux <fi as one can see 
in Figs. 9 and 10 showing the spectrum for different a's and gap angles 3 ; in 
the full-ring case we present a comparison with the homogeneous field having 
the same flux through V. Another difference from the homogeneous field case 
is that the discrete spectrum may be void for some values of the parameters 
which is again clear from the said figures; for the full ring the critical value 
is given by (|3.2|) . 
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3 The same is true in Fig. 8 where, however, the jump in the derivative is too small to 
be visible. 
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Figure captions 

Figure 1 The dependence of energy eigenvalues of a full ring on the 
radius R, a — 5. Each level is labelled by to, all levels except the 
ground state are twice degenerated. 

Figure 2 This figure shows a transition between the two regimes as the 
distance 9 between the loose ends of the ring increases. The full-ring- 
type spectrum (9 ~ 0) passes to the infinite-well-type around a ~ 29R. 
In the presented situation R = 2 and a = 10 is chosen. The dotted 
lines show the spectrum of a Dirichlet case. 

Figure 3 The contour plots of the (real) wave functions of the ground 
state and the 8th excited states of a broken ring are shown for a — 1, 
R = 10, and 9 = ir/3. The corresponding eigenenergies are E Q = 
-0.249 and E 8 = -0.00415, respectively. 

Figure 4 The influence of random choice of a on localization of the 
ground-state eigenfunction is presented. Each of the 5000 points rep- 
resents a particular (random) choice of the piecewise constant function 
a((p) that is formed by 10 constant pieces of the same angle length 
7r/5 and the value of a is uniformely distributed around ao = 1 with a 
given mean square Aa. The quantity (A^) is the (minimized) second 
moment of the probability distribution. The radius R = 5. 

Figure 5 The dependence of energy eigenvalues of a full ring on intensity 
of homogeneous magnetic field B, a — 1, R — 5. We plot only levels 
"setting on" the first three Landau levels. 

Figure 6 The dependence of energy eigenvalues of a broken ring on the 
gap angle, B = 0.2, a — 1, R — 5. Only first seven levels of each 
"band" are plotted. The dotted lines mark the Landau levels for the 
given B. 

Figure 7 The contour plots of the absolute value of wave functions of 
broken rings, B = 0.2, a — 1, R — 5. The upper subplot shows the 
ground state (E = -0.193) of a "almost full" ring 9 = vr/10 and the 
lower one an excited state (E = 0.552) of the second "band" with 

9 = 7T. 
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Figure 8 The dependence of energy eigenvalues of a full ring in the 
Aharonov-Bohm case on the magnetic flux 0, a = 1, R = 10. 

Figure 9 The dependence of energy eigenvalues of a full ring in the 
Aharonov-Bohm case on the magnetic flux 0, R = 1. The three val- 
ues of a = 0.5, 1, 1.5 are chosen so as to demonstrate both the non- 
smootheness of energy curves at integer values of and the possibility 
of a void spectrum for some intervals of 0. 

Figure 10 The comparison of dependence of eigenenergies on the mag- 
netic flux for a full 9 = and broken 9 = n/3 rings in the Aharonov- 
Bohm case, a = 0.5, R — 1. The dotted curve shows an analogical 
situation with the homogeneous magnetic field of the same flux. 
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